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The superconducting magnet system of a fusion reactor plays a vital role in plasma confinement, a process 
that can be disrupted by various operational factors. A critical parameter for evaluating the temperature margin 
of superconducting magnets during normal operation is the nuclear heating caused by D-T neutrons. This study 
investigates the impact of nuclear heating on a superconducting magnet system by employing an improved 
analysis method that combines neutronics and thermal hydraulics. 

In the magnet system, toroidal field (TF) magnets are positioned closest to the plasma and bear the highest 
nuclear-heat load, making them prime candidates for evaluating the influence of nuclear heating on stability. 
To enhance the modeling accuracy and facilitate design modifications, a parametric TF model that incorporates 
heterogeneity is established to expedite the optimization design process and enhance the accuracy of the compu- 
tations. A comparative analysis with a homogeneous TF model reveals that the heterogeneous model improves 
accuracy by over 12%. Considering factors such as heat load, magnetic-field strength, and cooling conditions, 
the cooling circuit facing the most severe conditions is selected to calculate the temperature of the superconduc- 
tor. This selection streamlines the workload associated with thermal-hydraulic analysis. This approach enables 
a more efficient and precise evaluation of the temperature margin of TF magnets. Moreover, it offers insights 
that can guide the optimization of both the structure and cooling strategy of superconducting magnet systems. 
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I. INTRODUCTION 


Superconducting magnets are widely used in energy, elec- 
tric power, transportation, medicine, and other fields [1]. 
NbTi and NbsSn are two of the most commonly used materi- 
als in superconducting tokamak fusion devices [2-5], which 
are cooled to temperatures of several K using liquid helium to 
achieve a superconducting state. However, during device op- 
eration, the superconducting magnet is likely to be disturbed 
by various factors such as mechanical vibrations, thermal dis- 
turbances, electromagnetic disturbances, flux jumps, neutron 
irradiation effects, and chemical corrosion. These types of 
disturbances may alter the superconducting-material compo- 
sition as well as its microscopic or macroscopic structure, 
resulting in the degradation of the superconducting perfor- 
mance of the conductor [6]. 

A superconducting magnet system is integral to a toka- 
mak fusion reactor and comprises toroidal field (TF) coils, 
poloidal field (PF) coils, and a central solenoid (CS) [7]. Dur- 
ing normal reactor operation, diverse thermal disturbances af- 
fect TF coils, with nuclear heating due to D-T neutrons play- 
ing a significant role [8, 9]. D-T fusion releases 4.52 MeV a 
particles and 14.06 MeV neutrons. While a particles are con- 
fined by magnetic fields, neutrons (electrically neutral) travel 
outward from the plasma area, interacting with local materials 
and leading to nuclear thermal deposition. Neutron—nucleus 
collisions within windings cause radiation damage, which de- 
grades the material properties when accumulated [10, 11]. 
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Hence, providing appropriate neutron shielding and cooling 
for superconducting magnet systems is crucial. Despite the 
substantial reduction in neutron flux via shielding materi- 
als such as blankets, divertors, shields, and vacuum vessels 
[12, 14], nuclear thermal deposition and neutron-irradiation 
effects remain substantial, necessitating a comprehensive sta- 
bility analysis. In thermal hydraulics, the cooling system 
must efficiently extract heat from the magnet system. Exces- 
sive heat deposition that surpasses the capacity of the cooling 
system leads to a localized rise in conductor temperature and 
eventual quenching [15]. Unlike PF and CS, TF coils primar- 
ily generate the main magnetic field; however, they are ex- 
posed to the highest neutron-flux levels [16]. Hence, rigorous 
neutronics and thermal-hydraulic analyses are imperative for 
TF coils to evaluate nuclear-heating disturbances and predict 
conductor temperature margins. 


Accurate assessment of the heat load on the TF coils re- 
quires the consideration of diverse factors such as nuclear 
heating, heat transfer from the case, cooling conditions, mag- 
netic fields, and strain. In the past century, projects such 
as ITER, JT-60SA, EU-DEMO, K-DEMO, and CFETR have 
conducted detailed calculations and analyses, frequently in- 
tegrating neutronics and thermal-hydraulic methods [17—26]. 
These methods have been optimized in several respects. JT- 
60SA’s nuclear-heating analysis [27] incorporated nuclear 
heating, coil-case heat transfer, and AC losses, noting a neg- 
ligible AC-loss contribution during normal operation. Simi- 
larly, ITER’s TF temperature-margin analysis [28] considered 
nuclear heating and heat transfer from the coil case, assess- 
ing the central and side pancakes under varying nuclear-heat 
load scenarios. DEMO [29, 30] and CFETR [9] employed 
two- and one-dimensional (2D/1D) thermal-hydraulic calcu- 
lations, considering the magnetic-field distribution, nuclear 
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Fig. 1. (Color online) Superconductor operating environment in fusion reactor. The red arrows represent nuclear heating on the TF case and 
WP, and the blue arrows represents the heat transferred from the case to the WP. 


heating, and coil-case heat transfer. 

In this study, an improved analysis method combining neu- 
tronics and thermal hydraulics is proposed to evaluate the 
influence of nuclear heating on the temperature of TF coils 
in a tokamak fusion reactor. Compared with the traditional 
method, the proposed method makes the following two im- 
provements: 

1) The heterogeneous parametric TF model accurately cap- 
tures the nuclear-heat distribution in cables, insulators, and 
other components. As fusion-reactor magnet-system de- 
signs evolve, this method adopts a parametric-modeling ap- 
proach, enabling the rapid adjustment of key parameters (such 
as cable diameter) for the quick generation of neutronics- 
calculation models, thereby saving time and enhancing ac- 
curacy. 

2) Compared with the conductor with the maximum mag- 
netic field, highest nuclear heating, and poor cooling con- 
ditions, only one or two cooling circuits are selected for 
thermal-hydraulics analysis. 

This method aims to precisely and efficiently assess the 
temperature margins during normal tokamak operations. Fig- 
ure | illustrates the TF-coil operating environment in a fusion 
reactor, and Fig. 2 outlines the proposed technical approach. 
The detailed steps are as follows. 

We define the analysis model, which encompasses the ref- 
erence CAD model, materials, and superconducting strand 
scaling parameters of the fusion reactor. A corresponding 
neutronics parametric heterogeneous model is created based 
on the CAD model, establishing cable, insulator, and jacket 
structures within the winding pack (WP) area. The impact of 
the improved method on the heat-load accuracy is presented 
in Sect. II] A. Neutron sources can be constructed using the 
plasma parameters (Sect. II B), and nuclear-data libraries such 
as ENDF/B [31], JEFF [32], and FENDL [33] are typically 
employed. The Nb3Sn/NbTi critical currents can be calcu- 
lated using the scaling-law equations. 


Monte Carlo-based neutronics calculations were used to 
assess the nuclear-heating distribution in TF coils caused by 
D-T neutrons and secondary photons. TF coils are typically 
poloidal segmented into several sectors to facilitate the statis- 
tical analysis of the nuclear-heating poloidal distribution. The 
conductor with the maximum nuclear heating was determined 
considering the radial, poloidal, and toroidal distributions on 
the WP. 

Thermal analysis software was used to calculate the heat 
transfer from the TF case to the WP using nuclear-heating 
data. The TF case was divided into a suitable number of 
poloidal segments, and a 2D finite-element model was em- 
ployed for the calculation. Different heat-transfer coefficients 
were applied to different boundaries. 

Neutronics calculations were combined with thermal- 
hydraulic analysis to calculate the heat load along the cooling 
circuit. The average nuclear-heating density for each turn was 
computed based on the material-volume proportions. The 
nuclear-heating density along the cooling circuits (such as 
HF-1) can be obtained by considering the cross-sectional area 
of each turn and the arrangement of the cooling circuits. 

One cooling circuit was selected for the thermal-hydraulic 
analysis based on the conductor with the maximum mag- 
netic field, highest nuclear heating, and poor cooling. A 
thermal-hydraulic analysis was performed considering the to- 
tal thermal-load, mechanical, cooling, and electromagnetic 
conditions, and the conductor-temperature evolution was cal- 
culated. 

This method offers a systematic and comprehensive evalu- 
ation of fusion neutron-induced TF-coil thermal stability, in- 
volving critical property calculations, neutronics modeling, 
nuclear-heating and heat-transfer analysis, conductor heat- 
load assessment, critical cooling-circuit analysis, and final 
temperature evaluation. The accuracy of the results was 
greatly improved with the elaborate neutronics model, and 
counting the radial and toroidal nuclear-heating distributions 
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Fig. 2. The technical approach of the method. The geometry and materials are input for the neutronics calculation, thermal analyses, and 1D 
thermal-hydraulics analysis. The heat-load calculation provides a heat source for stability analysis. 


was more convenient. Heterogeneous parametric neutronics 
models significantly boost the accuracy of the results and sim- 
plify the model modification process. This method facilitates 
a precise evaluation of the temperature margins for conduc- 
tors in the most challenging cooling circuit. Furthermore, 
it provides insights into optimizing the structure and cooling 
scheme of fusion-superconducting magnet systems. 


Sections II-IV demonstrate this method by analyzing TF- 
coil nuclear-heating stability in a tokamak fusion reactor. The 
major and minor radii of the reactor were 3.6 m and 1.1 m, re- 
spectively, with TF magnets generating a peak magnetic field 
of 13.5 T and central magnetic field of 6.5 T. The reactor 
was operated in the steady-state and transient modes. For 
the steady state, the fusion power was 40 MW with a 1000s 
plasma pulse. For the transient modes, the fusion power was 
180 MW with a 10s plasma pulse. The data were sourced 
from the corresponding magnet-design reports. 


I. ANALYSIS MODEL 
A. Conductor and strand parameters 


The TF-magnet system comprises 16 TF coils, each con- 
taining 126 turns. It carries a nominal current of 55kA and 
provides a 6.5 T magnetic field at the plasma center to con- 
fine the plasma. The WP is enclosed in a stainless-steel case. 
The geometry of the cable-in-conduit conductors (CICCs) is 
shown in Fig. 3(a). The WP can be divided into two sub-WPs: 
high-field (HF) and low-field (LF). The HF WP comprises six 
double pancakes (DPs), in which a high-Jc NbsSn strand is 
selected, whereas the LF adopts an ITER-like Nb3Sn strand. 
The structure of the conductor includes an insulator, jacket, 
cable, and helium channel, corresponding to gray, white, yel- 
low, and black, respectively, in Fig. 3(a). SS316LN was used 
for the TF case and jacket, and G10 was used as an insula- 
tor [34]. 

Figure 3(b) shows the horizontal profile of the cooling 
scheme for the HF WP. Six parallel cooling circuits are in- 
cluded, with two pancakes set as one cooling circuit. The 
length of each cooling circuit was approximately 80m. He- 


lium flowed in opposite poloidal directions in two adjacent 
cooling circuits. The primary parameters of the conductor 
and strands are listed in Table 1. 


The critical current in the Nb3Sn composites was calcu- 
lated according to the scaling law of [35]. The scaling pa- 
rameters of the Nb3Sn strand are listed in Table 2. 


B. Neutronics model 


Considering the symmetry of the fusion reactor, a sector 
spanning 22.5 degrees enclosed by reflective planes was used 
for the Monte Carlo neutronics simulations. The model en- 
compasses components such as the blanket, divertor, vacuum 
vessel, thermal shield, TF coils, PF coil, central solenoid, 
ports, and cryostat. 


The blanket, divertor, and vacuum vessel collectively con- 
tribute to neutron shielding, with the blanket assuming the 
primary role of shielding neutrons in the magnet system. As 
illustrated in Fig. 4, the TF coils were divided into 45 poloidal 
sectors arranged counterclockwise from the inboard leg. Sec- 
tors 12-23 are positioned at the rear of the divertor, whereas 
the remaining sectors are situated at the rear of the blanket. 
Sector 19 houses the helium inlet, and Sector 18 contains the 
outlet, as shown in the diagram. 

Owing to the intricate nature of the WP and the substan- 
tial number of torus surfaces, significant effort has been in- 
vested in neutron modeling using the constructive solid ge- 
ometry (CSG) method. The complexity of modeling in this 
area is typically addressed by employing a hierarchical struc- 
ture or homogeneous model. In line with the design of the 
TF coils, this study modeled a heterogeneous parametric WP 
model with a cable, stainless-steel jacket, and insulator struc- 
ture. 

However, the simulation faces limitations owing to the 
constrained capabilities of the Monte Carlo computing code 
MCNP [36, 37] when applied to a multitude of torus sur- 
faces. To maximally enhance model accuracy, an “approx- 
imate rectangular cross-section" was employed as a substi- 
tute for the “circular cross-section," determined based on 
material-volume proportions of the turns. 

Based on Fig. 5, both a homogeneous WP and a more in- 
tricate WP configuration were constructed to facilitate a com- 
parison for the precision enhancement. 

The probability distribution for sampling 14.1 MeV neu- 
trons is defined as follows: 


s(a) = [1 — (a/A)?]" ,0<a< A, (1) 


where the value of P is determined through physical mod- 
els and assumptions, typically ranging between 3 and 4. The 
parameter A corresponds to the minor plasma radius, and pa- 
rameter a defines the plasma contour. The representation of 
the magnetic surface in polar coordinates is expressed as fol- 
lows: 


R= Ro + a- cos(a + ô- sin a) 
+esh |1 — (a/A)?] (2) 
Z=El-a-sina 


Here, esh denotes the Shafranov factor, indicating the dis- 
placement from the magnetic-surface center to the low-field 
side. Ro is the major plasma radius. Æl and 6 are the elonga- 
tion ratio and triangularity of the plasma, respectively [38]. 

During this phase, the principal parameters of the conduc- 
tor and strand are initially established, and the critical current 
in the Nb3Sn composites is computed using scaling-law equa- 
tions. Subsequently, a heterogeneous parametric neutronics 
model and neutron source were constructed. These models 
were based on the reference CAD model and plasma param- 
eters, respectively. 


I. HEAT LOAD 


The cumulative heat load on the conductors comprised nu- 
clear deposition (from both neutrons and photons emitted by 
D-T reactions) and heat transferred from the case to the WP. 
A higher nuclear deposition on the central turn than on the 
corner turn can be reasonably anticipated. However, the sit- 
uation becomes more intricate when the heat transfer from 
the case is considered. The influence of displacement and 
gaps on the heat transfer in the TF coils under Lorentz forces 
during normal operation must be considered [16]. Previous 
research [39] has highlighted that corner turns experience a 
notably greater heat-transfer load than central turns. Con- 
sequently, a precise statistical analysis of both toroidal and 
radial heat load profiles is essential. 

By comparing the magnitudes of the nuclear heat and the 
heat transferred across various turns, cooling circuit capable 
of effectively managing the maximum heat load may be iden- 
tified. 


A. Nuclear heating 


Nuclear heating was evaluated using the Monte Carlo 
particle-transport code MCNP v5.1.6 with the FENDL 3.1d 
nuclear-data library. Assuming a fusion power of 40 MW, 
the corresponding D-T neutron-emission rate was set as 
1.42 x 101° s71. A bi-dimensional (radial and poloidal) dis- 
crete neutron source was employed to characterize the spa- 
tial distribution of D-T neutrons. By employing the neu- 
tronics model above, the nuclear-heating distributions within 
the TF coils were calculated based on 10’ neutron histo- 
ries. To enhance the statistical efficiency of the Monte Carlo 
simulations for shielding calculations, an importance-based 
variance-reduction technique was implemented to ensure sat- 
isfactory results and statistics, which generally remained be- 
low 0.5%. 

Table 3 outlines the nuclear-power deposition in the TF 
coils. The overall TF-coil deposition amounts to 17.86 kW, 
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Fig. 3. (Color online) (a) A half-horizontal profile of the TF coil at the inboard leg half. The numbering of the pancakes and the turns adopted 


in this work are shown. (b) Cooling scheme of high field. 


Table 1. Main conductor and strand parameters 


HF WP conductor LF WP conductor 
Strand type OST-E2013 ITER-like 
Strand diam. (mm) 0.82 0.82 
No. of SC 720 450 
No. of Cu 420 300 
Cable layout (2SC+1Cu)x3x4x5x6+lcore ((2SC+1Cu) x3x3x5+1core)x5 
Cu core layout 3x4x5 3x5 
Cable pattern STP STP 
Jacket material 316LN 316LN 
Void fraction ~30% ~30% 
Cable diam.(mm) 34.3 27.6 
Axial strain -0.6% -0.6% 


Table 2. Scaling parameters of Nb3Sn strand 


Parameters High-Je _ ITER-like 
Nb3Sn  Nb3Sn 
Ga 49.23 49.00 
Caz 7.70 0.300 
Sis 0.32% 0.312% 
Emax —0.07% 0.059% 
Bčzwmax(T) 31.68 33.24 
Tht) 1653 16.34 
C(AT) 58158 21700 
p 0.60 0.593 
q 2.17 2.156 


with the case contributing 51.7% and WP providing 48.3% 
of this total. Additionally, a comparison was made between 
the nuclear-power deposition at the rear of the blanket and di- 


vertor with the aim of evaluating the neutron-shielding effec- 
tiveness. Notably, the total nuclear heat of TF coils obtained 
through the homogenized model reached 17.46 kW, exhibit- 
ing a relative error of 2.2% compared with the detailed hetero- 
geneous model. Considering the statistical-count error (below 
0.5%), the influence of the model accuracy on the comprehen- 
sive nuclear heat cannot be dismissed. 


To further assess the precision enhancement of the hetero- 
geneous model in terms of TF nuclear heating, a comparison 
of two typical cross-sections is shown in Table 4: one at the 
inboard middle plane (IB) and the other at the outboard mid- 
dle plane (OB). The differences in the nuclear-heating densi- 
ties of the cable, jacket, and insulator in the IB were —16.2%, 
7.3%, and 2.7%, respectively. Significant differences were 
observed in the OB. This is because in the heterogeneous 
model, only neutrons passing through the case, jacket, and 
insulation interacted with the cable material, leading to the 


Table 3. Nuclear power deposition of TF coils (kW) 


Poloidal sector WP-cable WP-jacket WP-insulator Case Total 
BLK rear 2.25, 2.38 0.77 4.93 10.34 
DIV rear 1.28 1.41 0.52 4.31 7.52 

Total 3.53 3.79 1.30 9.24 17.86 
Percent 19.76% 21.22% 7.28% 51.74% 100% 


BLK ear are sectors of TF coils behind the blanket module, DIV rear are sectors of TF coils behind the divertor. 


Table 4. Nuclear heating densities at inboard middle plane (IB) and outboard middle plane (W / cm?) 


Model IB-cable IB-jacket IB-insulator OB-cable OB-jacket OB-insulator 
Hetero. model 7.57 x 107 9.84 x 107 3.56 x 107 1.63 x 107 2.09 x 107 8.61 x 10-° 
Homo. model 8.80 x 107+ 9.12 x 107+ 3.46 x 107+ 1.83 x 1074 1.95 x 1074 8.32 x 107° 

Difference —16.20% 7.30% 2.70% —12.50% 6.70% 3.40% 


Difference = (Value of heterogeneous model- Value of homogeneous model)/Value of heterogeneous model x 100%. 
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Fig. 4. Poloidal sectors of TF coils. 


Fig. 5. (Color online) Horizontal cross section of the inboard neu- 
tronics model. 


observed variations. In the thermal-stability analysis of the 
TF-magnet system, the precision of the cable in the hetero- 
geneous model improved by more than 12%. This improve- 
ment is crucial for subsequent heat-loading and temperature- 


margin assessments in the thermal-hydraulic nuclear context. 


Figure 6 shows the polar distribution of the WP nuclear- 
heating density within HF conductor layer 1. The average 
nuclear-heating densities of the cables, jackets, and insula- 
tors were 3.35 x 1074 W/cm, 4.34 x 1074 W/cm, and 
1.63 x 1074 W/cm, respectively. For the nuclear-heating 
density, the insulator displayed the lowest value, the cable 
registered approximately 2.06 times the value of the insulator, 
and the jacket demonstrated a nuclear-heating density approx- 
imately 2.67 times greater than that of the insulator. The high- 
est value was observed on the upper side of the IB (Sector 5). 
Furthermore, the conductor (Sector 16—18) at the rear of the 
divertor exhibited an elevated nuclear-heating density, sug- 
gesting its potential for divertor-shielding enhancement. Neu- 
tron fluxes in proximity to the middle plane and port tended to 
be high, accompanied by corresponding nuclear-heating den- 
sities on the rear side of TF coils. 


Figure 7 illustrates the exponential decay of the nuclear- 
heating density along the radial direction for the WP in Sec- 
tors 5 (IB) and 30 (OB). Comparing the rear layer in the HF 
with the front layer at the IB, the value decreased by 3.56 
times in the HF and by 16.26 times in the LF. Similar com- 
parisons with the OB front layer showed a 3.14 times decrease 
in the HF and only a 4.99 times decrease in the LF. Notably, 
the nuclear-heating density in the HF surpassed that in the LF 
by an order of magnitude, highlighting the preference for the 
HF in subsequent analyses. 


The statistical assessment of the nuclear-heating toroidal 
profile in the HF is important for selecting the optimal cool- 
ing circuit for thermal-hydraulic analysis. As shown in Fig. 8, 
a comparison between Sectors 1 and 5 revealed that, in Sec- 
tor 1, the cable, jacket, and insulator values in HF-6 (cen- 
tral pancakes) were 13%-16% higher than those in HF-1 and 
HF-2 (side pancake). A similar trend was observed in Sector 
5, with an increase of approximately 5%—-8%. As expected, 
the central turn exhibited a higher nuclear deposition than the 
corner turn; however, the discrepancy was not notable. 


The TF case comprised U-shaped enclosures (AU and BU) 
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Fig. 6. Nuclear-heating density poloidal profile in first layer of WP. 
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Fig. 7. Nuclear-heating density radial profile of the WP in Sector 5 (IB) and Sector 30 (OB). 


Table 5. Heat flux of several turns in Sector 5 


Turn HF 1-1 HF 1-2 HF 1-3 HF 1-4 
Heat flux (W/cm) 8.063 x 1078 3.663 x 107% 3.223 x 1073 3.916 x 107° 
Turn HF 2-1 HF 2-2 HF 2-3 HF 2-4 
Heat flux (W/cm) 3.289 x 1073 1.287 x 1073 1.287 x 1073 1.045 x 1073 
Turn HF 6-1 HF 6-2 HF 6-3 HF 6-4 
Heat flux (W/cm) 3.256 x 107% 1.287 x 1073 1.287 x 1073 1.023 x 107% 
B. Heat transferred from case to coil cooling channels, whereas AP and BP had eight cooling chan- 


nels each. The mass-flow rate for each cooling channel was 


and cover plates (AP and BP). Both AU and BU featured 22 


2.7 g/s, with an inlet pressure of 6 bar and inlet temperature 
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Fig. 8. Nuclear heating toroidal profile of the WP in Sector 1 and Sector 5. 


of 4.2K. During operation with a pulse duration of 1000 s, 
a Steady state for heat transfer between the case and conduc- 
tor was achieved. However, in the case of a pulse duration 
of 10 s, heat transfer was disregarded. Ansys 19.0 was uti- 
lized to conduct a 2D heat-transfer analysis, relying on the 
nuclear-heating results obtained from MCNP for the case at 
40 MW. Owing to the inward compression resulting from the 
centering force on the WP, a millimeter-sized gap was formed 
between the inner side of the WP and the case [16]. Conse- 
quently, the heat transfer between the case and plasma-side 
of the WP became limited. An ideal heat transfer between 
the case and other sides of the WP is assumed considering 
the isotropic material properties. The heat-transfer coefficient 
was calculated using the following equation: 


hra = 0.027 x Re®® x Pr®* x d/ Dh, (3) 


where Re and Pr represent the Reynolds and Prandtl num- 
bers, respectively. A denotes the thermal conductivity, Dy is 
the hydraulic diameter, and ha represents the heat-transfer 
coefficient. 

The finite-element model of the TF case for the heat- 
transfer calculation was discretized into 45 2D sectors in the 
poloidal direction. As indicated in Table 5, the heat flux trans- 
ferred from the TF case to the conductors was calculated. No- 
tably, corner turns exhibited a significantly higher heat flux 
than central turns. In the toroidal direction of Sector 5, the 
heat-flux value in HF 1-1 was approximately 2.5 times greater 
than those in HF 2-1 and HF 6-1. In the radial direction, the 
heat-flux values in HF 1-1, HF 2-1, and HF 6-1 surpassed 
those in the turns of the second layer by a substantial margin. 

These outcomes highlight the substantial variation in heat 
transfer from the case to the coil across different turns, and 


these values constituted a significant portion of the overall 
heat load. 


C. Total heat load on conductor 


In the subsequent phase, the cumulative heat load absorbed 
by the conductors along the cooling circuits is quantified. 
This process can be divided into two steps: first, the average 
nuclear-heating flux for each turn is calculated based on the 
material proportion and cross-sectional area of the turn; sec- 
ond, the nuclear heat on HF-1 and HF-6 is determined based 
on the arrangement of the helium cooling circuits. The cal- 
culation of the average nuclear-heating flux is given by Eq. 
(4), where q and V symbolize nuclear-heating density and 
the volume proportion of the cable in each turn, respectively. 
A represents the cross-sectional area of the turn. 


N Hux = (deable i Veable T jacket * Vacket T 


4 
z Vinsulator ) s Atum ( ) 


insulator 


In this scenario, the helium inlet and outlet were strategi- 
cally positioned in Sectors 19 and 18, respectively. Conse- 
quently, the heat load along the cooling circuits within each 
of the four layers (each approximately 420 m) was initiated 
at the heat flux originating from Sector 19 and concluded in 
Sector 18. 

Figure 9 shows the heat load on the conductor along the 
cooling circuits. HF-6 HT (heat transferred from the case 
to HF-6) was approximately half that of HF-6 NH (nuclear 
heating on HF-6). While the nuclear heating on the side pan- 
cake (HF-1 NH) was marginally lower than that on the central 
pancake (HF-6 NH), the heat load transferred from the case 
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Fig. 9. (Color online) Nuclear heating (NH), heat transferred (HT) and total heat load (Total) along HF-1 and HF-6 cooling circuits. 


(HF-1 HT) significantly surpassed that of the central pancake 
(HF-6 HT). Thus, HF-1 experienced the highest heat load, 
followed by HF-6. 

Notably, the analysis did not consider the impact of AC 
losses on the heat load. Owing to the unique cooling scheme 
and variations in the heat-transfer coefficient and WP struc- 
ture, the heat load of the central and side pancakes may dif- 
fer. Thus, a comprehensive analysis of the maximum heat 
load that considers both the distribution of nuclear heating 
and conduction heat of the coil box is imperative. 

Furthermore, the conductor heat-loading method can be ap- 
plied to alternative arrangements of helium cooling circuits. 
Of particular significance are the TF sector numbers, where 
the helium inlet and outlet are located, and the poloidal direc- 
tion of the helium flow. 


IV. THERMAL-HYDRAULIC ANALYSIS 
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Fig. 10. (Color online) Selection of helium cooling circuits. 


In this phase, the cooling circuits subjected to the peak 
magnetic field and heat load were examined, and the cooling 
circuit with the most challenging combined conditions was 
selected for 1D thermal-hydraulic analysis. The simulation 


results were obtained using Gandalf [40], a widely used tool 
for magnet-stability analysis. The input parameters for the 
thermal-hydraulic analysis are listed in Table 6. 


The choice of the cooling circuit is depicted in Fig. 10. 
HF1-1 and HF12-1 experienced a peak heat load; HF6-1 and 
HF7-1 encountered a peak magnetic field and high heat load; 
and HF7-1 represented the turn facing the peak magnetic 
field under poor cooling conditions. This is because of the 
clockwise helium flow in HF6-1 and counterclockwise flow 
in HF7-1. This caused hydraulic length from Sector 19 (he- 
lium inlet) to Sector 5 (peak magnetic field and nuclear heat- 
ing) in HF7-1 to be longer than that in HF6-1. Considering 
the magnetic field, heat load, and cooling conditions, HF7 
was chosen for the thermal-hydraulic analysis. It comprised 
four turns, forming an 80-m path. 


The thermal-hydraulic analysis model comprised four 
parts: the strands, conduit, bundle helium, and hole helium. 
The convective heat-transfer calculations were based on a pre- 
vious study [41]. Adiabatic boundary conditions were applied 
to the ends of the cables. 


Three aspects of TF operating conditions were considered: 


(1) Mechanical conditions: The influence of the Lorentz 
force on the superconducting magnet operation necessitates 
the consideration of strain. The axial strain of the supercon- 
ducting strand was assumed to be —0.6% under these calcu- 
lations. 


(2) Cooling conditions: Assumed cooling conditions en- 
compassed an inlet temperature of 4.2 K, pressure of 6.0 bar, 
and mass-flow rate of 8.0 g/s. The hydraulic length of the HF 
was approximately 80 m. 


(3) Electromagnetic conditions: The conductor carried a 
steady current of 51.7 kA, resulting in a peak magnetic field 
of 13.5 T on the WP. 
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Table 6. Input parameters of the thermal-hydraulic analysis 


Parameters Value 

NbsSn cross section (m7) 1.90 x 1077 
Copper cross section (m°) 4.12 x 1074 
Jacket cross section (m°) 6.36 x 107+ 
Insulation cross section (m?) 2.29 x 1074 
Helium in Bundle cross section (m?) 2.94 x 107+ 
Helium in Hole cross section m’? 2.83 x 107° 
Bundle-region hydraulic diameter (m) 4.74 x 1074 
Hole-region hydraulic diameter (m) 6.00 x 107° 

Wetted Perimeter between strands and the bundle helium (m) 2.45 
Wetted Perimeter between conductor strands and the jacket (m) 1.08 x 107? 
Wetted Perimeter between jacket and the bundle helium (m) 9.70 x 107? 
Wetted Perimeter between bundle and the hole helium (m) 1.88 x 107? 

Percentage perforation 0.1 


Temperature (K) 


Distance from helium inlet (m) 


Fig. 11. Temperature evolution of the superconducting strand in steady-state mode. 
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Fig. 12. Conductor temperature and pressure history at the helium outlet. 


A. Steady-state mode 


In the steady-state mode, both nuclear heating and heat 
transferred from the case to the coil were considered. Fig- 


ure 11 illustrates the temperature evolution of the supercon- 
ducting strand. The temperature history at the helium-flow 
outlet and pressure decrease are shown in Fig. 12. The 
conductor-load temperature sequentially increased with dis- 
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Fig. 13. (Color online) Temperature evolution of the superconducting strand in transient mode. 


tance from the helium inlet owing to the heat load, and the 
average temperature continued to increase over time. By the 
end of the 1000s plasma pulse, the conductor temperature 
at the end of the helium-flow channel reached 5.502 K, indi- 
cating an increase of 1.3K from the initial value of 4.2 K. 
The most substantial temperature increase occurred in the 
first turn (0 m—20 m distance from the helium inlet) of HF-7, 
contributing approximately 45%-50% of the total tempera- 
ture increase. A milder increase was observed in the last turn 
(60 m—80 m distance from the helium inlet). The conductor 
temperature at the helium outlet reached its peak at 550s, 
started to decrease after 1000s, and reached its minimum 
value after1550s. The outlet pressure was 5.536 bar. The ini- 
tial temperature margin under the aforementioned peak mag- 
netic field and strain conditions was 3.0 K. After subtracting 
other heat loads apart from nuclear heating (such as hystere- 
sis loss, coupling currents, thermal radiation, and conduction, 
which amount to approximately 0.45 K), a remaining temper- 
ature margin of 1.25 K was achieved, satisfying the design 
requirement of > 1.0 K. 


B. Transient mode 


In the transient mode, the nuclear heat directly absorbed by 
the conductor remained constant throughout the 10-s pulse. 
However, the heat transferred from the case was not included 
in the analysis. Unlike the extended plasma-pulse mode, 
where the conductor and coolant can approach a thermal- 
equilibrium state, the 10-splasma pulse prevents them from 
achieving such an equilibrium. Consequently, in this mode, 
the total heat load can be approximated as 4.5 times the nu- 
clear heating load in the steady-state mode. 

Figure 13 shows the temperature evolution of the super- 
conducting strand over a 10-s pulse. In general, the conduc- 
tor temperature followed an increasing trend that aligned with 
the nuclear-heating load. The first turn of HF-7, which expe- 


rienced the greatest nuclear heating, exhibited the most sig- 
nificant temperature variation. At the end of the pulse, a peak 
temperature of 4.583 K was reached, indicating an increase in 
0.38 K over the inlet temperature. After deducting the other 
thermal loads, the remaining temperature margin was approx- 
imately 2.17 K, which satisfied the design requirements. The 
outlet pressure was maintained at 5.536 bar. Compared with 
the steady-state mode, the transient mode afforded a larger 
temperature margin. 


V. DISCUSSION AND SUMMARY 


The proposed method effectively combines neutronics and 
thermal-hydraulic analyses, resulting in a comprehensive as- 
sessment of the effect of the heat load, magnetic field, and 
cooling conditions on TF coils. Compared with traditional 
methodologies, this improved method introduced two crucial 
advancements in neutron modeling and thermal-hydraulic 
analysis. 

To demonstrate the application of this improved approach, 
an under-designed fusion reactor was used. During the 
neutronics-modeling phase, a parametric modeling technique 
was employed to construct a heterogeneous TF model. This 
not only reduced modeling time for model changes in the 
design-optimization process, but also enhanced the counting 
accuracy by more than 12%. Before carrying out a conductor- 
temperature analysis, the proposed method considers the cu- 
mulative effects of total heat load, magnetic field, and cooling 
conditions. Subsequently, one or two cooling circuits facing 
the most severe conditions were identified for the subsequent 
thermal-hydraulic calculations. This streamlined the work- 
load for subsequent analyses. 

Building on the aforementioned method, a more precise 
analysis of the thermal impact on TF temperature under the 
standard operational conditions of a fusion reactor can be 
achieved, thereby expediting the design process. Further- 


more, this method offers pivotal insights for optimizing the 
structure and cooling arrangement of fusion superconducting 
magnet systems. 
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